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Abstract 

We examine the dynamics of a self-gravitating domain wall using the 
model as a specific example. We find that the Nambu motion of the wall is 
quite generic and dominates the wall motion even in the presence of gravity. 
We calculate the corrections to this leading order motion, and estimate the 
effect of the inclusion of gravity on the dynamics of the wall. We then treat 
the case of a spherical gravitating thick wall as a particular example, solving 
the field equations and calculating the corrections to the Nambu motion an- 
alytically for this specific case. We find that the presence of gravity retards 
collapse in this case. 
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I. INTRODUCTION 



Over the past few decades, topological defects have become a familiar class of objects 
in many areas of physics. In cosmology, defects are believed to arise generically during 
phase transitions in the early Universe, and have been notably invoked to account for the 
anisotropics which seeded cosmological structures Despite the discouraging discovery 
that the power spectrum of global strings does not agree with observations of the COBE 
satellite 0, recent claims of a nongaussian component in the microwave background 
(see also [^]), and the improved agreement of the spectrum predictions for models of cosmic 
strings with a cosmological constant |^ (whose best fit coincides with the value of A deter- 
mined from type la supernovae P]) suggest that it may be still too early to discard them as 
the source of cosmic structure. 

Domain walls are defects that arise when the phase transition occurs by the breakdown 
of a discrete symmetry. They correspond to solitons in 1 + 1 dimensions which are extended 
in two spatial dimensions to form a wall structure. Because static wall solutions depend only 
on one coordinate (the distance from the wall's core), they can often be found analytically 
in the absence of gravity, and perturbatively analytically in its presence. In a cosmological 
context, it was soon realised [0] that the existence of domain walls with ri>l MeV must be 
ruled out, because a network of such defects would rapidly evolve to dominate the energy of 
the Universe. Nevertheless, domain walls remain intrinsically interesting objects to study, 
for instance for their properties as hyp ersurf aces; in cosmology, domain walls have been 
proposed as a realization of our universe in higher dimensions P], and are currently being 
explored as a possible resolution of the hierarchy problem [§. 

When considering 'defects', there are two main aspects to understand: their gravita- 
tional (or other particle) interactions and their dynamics. The gravity of domain walls is 
an interesting and rather more subtle topic than it might seem at first sight. Indeed, unlike 
all other defects (with the exception of global strings JlOl)) the wall's metric is not in gen- 
eral static |Tl|JT^ , but admits a de Sitter-like expansion in its plane. Moreover, observers 
experience a repulsion from the wall, and there is a cosmological horizon at a finite proper 
distance from the defect's core. This horizon is a consequence of the choice of coordinates, 
and in a different set of coordinates |jl3[ the wall has the appearance of a bubble which con- 
tracts in from infinite radius to some minimum radius, then re-expands, undergoing uniform 
acceleration from the origin. The 'horizon' is then simply the lightcone of the origin in these 
coordinates, and is somewhat similar to the horizon of Rindler spacetime. These results 
were originally obtained for infinitesimally thin walls, using the hypersurface formalism de- 
veloped by Israel [|14|, but can be shown to be robust as an approximate description of a 
thick domain wall by a perturbative expansion in the thickness of the wall []T3|,|TC[ , or within 
the context of a fully nonlinear treatment of a scalar field coupled to gravity fl^ . 

The crucial physical difference, then, of the self-gravitating domain wall spacetime, is the 
presence of the cosmological horizon, which introduces a second length scale into the system. 
Ordinarily, a defect possesses one length scale, its thickness w, however, the distance to the 
event horizon of the domain wall gives another length scale, Wh, which can be compared to 
w. These lengths are given in terms of the coupling constants of the theory and, as taking 
a thin wall limit turns out to be a very artificial construction in terms of these underlying 
parameters, it becomes pertinent to examine both the gravity and dynamics for a thick 
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domain wall. 

The dynamics of topological defects are typically extremely nonlinear, and their study 
is usually carried out in the so-called Nambu approximation, where the full-field theory 
action of the model in which the defects arise is replaced by an approximation based on the 
degrees of freedom of the defect's core. At first sight, this looks like a perfectly reasonable 
approximation, since — even for finite-sized defects, such as string loops and wall bubbles — 
the thickness of the defect is typically many orders of magnitude smaller than its size. It 
must be noted however, that the Nambu action is only a leading order approximation of the 
real action, and is obtained in the limit a — > 0, where a oc Kw is proportional to the defect's 



typical curvature and thickness |18-|21|. However, it is often the case that the small-scale 



structure of defects (when a can be significant) is of particular importance when considering 
their impact. It seems that quite generically these are the points at which the defects lose 
most of their energy (by Higgs, gravitational or possibly other types of radiation, depending 
on the type of defect) and the properties of the small-scale therefore have a direct impact 
on the defects' lifespan, which in turn dictates the possible cosmological implications they 
can have. 

The first attempts to derive effective actions for walls p!9|-pl| naturally neglected gravity, 
in the sense that neither curvature of the background spacetime nor the self-gravity of the 
wall were considered. The effect of background spacetime curvature was considered (using 



the effective action method discussed below) in [^, and the self-gravity of the wall for a 
very special trajectory was considered in [^]; however, the motion of a fully self-gravitating 
thick domain wall (where the curvature of spacetime is that induced by the wall) has not to 
date been considered. In this paper, we address this problem, namely, using the A$'^ field 
theory, we examine the dynamics of a thick, self-gravitating kink domain wall solution. 

There are essentially two methods which have been used to obtain the effective motion 
of a thick defect, both of which involve in some way an expansion of the fields around 
a well-known solution (such as the hyperbolic tangent kink for a A$'^ wall). The first 
approach consists of replacing the solution in the action and integrating out perpendicular 
to the defect, which yields an effective action based on the defect's core. This method 
has been employed for instance in p8|j20| -p2|. The second approach consists of examining 



the field equations perturbatively in a relevant parameter (or parameters) such as a, with 
the equations of motion to a particular order arising as an integrability condition (see for 
instance ||20|). Here we use this latter method, adapted to take the wall's self-gravity into 
account. The reason for this choice is in fact related to the inclusion of gravity in the 
problem; the motion of the core of the defect interacts with massless degrees of freedom in 
the bulk (the graviton) and hence a correct application of the effective action method is less 
transparent. 

The layout of the paper is as follows. In the next section we present the model, then 
briefly introduce the Gauss-Codazzi formalism and derive the Einstein equations in the 
corresponding "3 + 1" notation. We end this section by reviewing our method in the case 
of a flat background spacetime. In section |ITT| we solve the fleld equations for the case a > e 
(where e characterizes the gravitational interaction of the scalar fleld). In section |^ we 
discuss the particular case of a collapsing spherical wall, and we conclude in the last section. 
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II. THE GAUSS CODAZZI FORMALISM 



Our starting point is the usual Goldstone matter Lagrangian, 

/: = (v„$)'-f/($), (1) 

where $ is a real Higgs scalar field and U{^) is a symmetry-breaking potential which we 
take to be 

[/($) = A ($2 - r]^y = Ar/V (^^^ . (2) 

This model admits domain wall solutions, where the Higgs field tends to different vacua 
$ = ±1] at, say, u = ±oo for a fiat wall. This implies the existence of a surface for which 
$ = 0, and this surface defines the defect's core. 

As in [|l^ , we scale out the dimensionful parameters from the Lagrangian by defining 

X = ^/rj e = S-KGrf] (3) 

X now tends to ±1 at the vacua. The scalar and Einstein equations are 

UX + ^X(X^-1) =0 (4a) 



T^ab = e 



2X aX b gab (x"^ — 1 



(4b) 



where TZab is the spacetime Ricci tensor, and w is proportional to the inverse mass of the 
Higgs after the symmetry breaking, and therefore characterizes the width of topological 
defects within the theory. Note that according to (^) e characterizes the gravitational 
interaction of the Higgs field. 

As mentioned in the introduction, our method for finding the effective dynamics of 
domain walls requires the expansion of the quantities appearing in the full equations of 
motion in powers of some small parameter. This can be achieved by splitting these quantities 
in their components parallel and perpendicular to the wall's worldvolume S, with the help of 
the Gauss-Codazzi formalism [^. Before we start, a remark on our notation: although we 
shall generally use lowercase Latin indices a,b, . . ., we may emphasize the parallel character 
of some index by using uppercase Latin letters A,B, . . . For instance, the coordinates parallel 
to the defect will be called a^. The proper distance from the wall will be denoted by u. 

A domain wall's core (defined by the location of X = in the above model) is a three- 
dimensional surface in four-dimensional spacetime, and consequently it admits a (spacelike) 
unit normal field denoted by ra". This normal field can be regularly extended off the world- 
volume by imposing n°''W = 0, so that each surface of constant u has a normal field n", 
a first fundamental tensor hab and a second fundamental tensor Kab, the latter two being 
defined by 

hab = gab + riaUb Kab = hf" aV cTib- (5) 

hab is the projection tensor onto the worldvolume (with our choice of coordinates, its parallel 
components are equal to the intrinsic metric's, and its perpendicular components vanish). 
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Kab is the extrinsic curvature of S, which hes tangentially to the worldsheet and describes 
how the wall curves away from a hyperplane in spacetime. The metric of the wall spacetime 
is therefore written in this coordinate system as 

(is^ = hABdcT^da^ — du^. (6) 

The Gauss and Codazzi integrability conditions for the hypersurface generated by the 
wall's core are, 

Rab = n^dKht + - Kkab + nacbdU'n'' (7) 



b,kl - Dak = RdeU^ht 



(8) 



where Da = h'^V c- Using these integrability conditions, the equations of motion (0) can be 
rewritten 



ab 



du 
du 



2KacKl - KKab 



R. 



ab 



dX 

-k— + DaD-X 
ou 



dK 

du 



^ab 

dX 



2X'2 + {X^-l] 



D,K^„ - D„K = 2e—DnX 



du 

2 

^ab 



R = -k^ + 2e 



X'2 - (X^ - 1)^ + DaXD^X 



(9a) 

2eDaXDbX (9b) 
(9c) 

(9d) 

(9e) 
(9f) 



Here we have (without loss of generality) set w = 1, which amounts to using wall rather than 
Planck units. The mean curvature K is the trace of Kab and Rab is the worldvolume Ricci 
tensor. Note that the equations for Kab, K and the derivatives of K in the wall correspond 
respectively to the "Ai?," "mm" and non-diagonal Einstein equations. The final equation 
results from the trace of (0) and turns out to be related to the integrability condition which 
gives the wall equation of motion. 

We now need to identify the two key parameters in these equations, as well as the de- 
pendence of the variables on these parameters. Clearly the gravitational parameter appears 
explicitly in the equations, however, the parameter a characterising the motion of the wall 
is only implicit in the equations. The first step to identifying this parameter is to quantify 
what one means by 'motion' of the wall; this is encoded in the components of the extrinsic 
curvature on the wall core itself, since this tells us that the wall is curved in the ambient 
spacetime. We therefore set 



a 



\Kt{u 



0)|. 



(10) 



Note that this is not a fundamental parameter of the theory, in that it is not given in terms 
of any coupling constants or masses, but simply represents the physical motion of the wall 
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and enables the effect of that motion to be correctly considered. The extrinsic curvature Kab^ 
therefore, has two main contributions: that of the motion of the wall core and that of the 
gravity of the wall. This can be estimated by considering the case of the plane-symmetric 
wall [|T^, for which a = 0, and 

i^,, = A'Adiag(l,-e2^*,-e2^*), (11) 

where A{u) = 1 + e + O(e^), c = 2e/3 + O(e^). Clearly then, the components of 

are 0(e). 

To summarise: in order to describe the motion of a domain wall in a curved spacetime we 
naturally have two parameters; one characterising the motion of the wall itself, a, and one 
the curvature of the ambient spacetime, e. The basic procedure for determining the equation 
of motion of the wall is to solve (|^) order by order in these parameters, investigating any 
constraints arising on the extrinsic curvature at each step. It is perhaps worthwhile briefly 
reviewing this process for e = 0, since the methodology is very similar when gravity is 
included. 

We begin by rescaling the extrinsic curvature and parallel coordinates via 

kab = aKab (12a) 
= x^/a. (12b) 

Since we are in flat space the Gauss identity (|^) simplifies to, 

Rab = KacKl - KKab, (13) 

hence 

Kab = Kab\o + auK^\oKbc\o (14) 

(where "|o" indicates that a quantity is evaluated at the wall core) is actually an implicit 
exact solution to the i^-equation. To order a, we see that 

K = K\o- auKl,\o, (15) 

and we can examine the X-equation ( P^ ) by setting X = X + aX-^, where Xq = tanhw, 
finding 

VX, = X'; - 2X,(3X2 - 1) = -K\oX'^. (16) 

Ordinarily, we might expect to be able to write the solution in terms of the basis of 
eigenfunctions of the operator V, however, we cannot do this directly, since X^ is in fact 
the zero mode of V. We are therefore forced to either deduce that K\q = 0, or, we can take 



the approach of reference |]2T| and remove the requirement that the equations of motion be 
regular at the wall. Since we are looking for freely moving wall trajectories, we will take the 
former approach, which can be summed up as an 'integrability requirement': multiplying 
both sides of (|l^) with Xq, and integrating over IR implies 

/oo 
Ko\o{Xl,rdu. (17) 
-oo 
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In order for Xi to have the appropriate asymptotic behavior for large u it follows that 



Ko\o = X, =0). (18) 

Equation (p!8| ) is of course the Nambu equation K\q = 0, which simply means that the core 
of the defect, situated at m = 0, follows the Nambu dynamics to zeroth order. 
We can then repeat this process, expanding order by order to get 

hab = Kblo + '^(yuKablo + a'^u^KacloK^lo (19a) 
K = K\o- auKl,\o + a\^K^\oK%K:\o (19b) 
X = Xo + aXi + a^X2, (19c) 

where K\q = O(a^), and 

ru ru 

X2 = sech^M / cosh u / {uK1j)^q) sech u. (20) 

JO "'-00 

This is sufficient to obtain the leading corrections to the Nambu action via the integrability 
constraints to third order. To third order (pH) gives 



hence 



[X',X', - X,X'^ = -K2X',' = -{K2\o + u'KS\oK\oK:\o)X',', (21) 



K2\o = -^^^K^\oK%K:\o = -(:-- 1 ) K^\oK%K^^\o, (22) 
/o(oo) \6 J 

where /„(«) = Jq duu'^ sech.'^u. 

Note that this process of using the integrability condition to derive a constraint on the 
extrinsic curvature uses the X-equation to order 0(a""'"^) for a constraint on K to order 
0(a"), and that this constraint only involves the even part of K, since any odd parts 
integrate to zero. Moreover, the i^-equation (^) shows that the even part of K to order 
0(a") depends on the odd part of K'^f^ to order 0(a"~^). Keeping this observation in mind 
prevents the unnecessary calculation of corrections to the geometric parameters. 



III. THE MOTION OF A WALL WITH GRAVITY 

In order to include gravity, we will make the initial assumption that gravity is subdom- 
inant to the motion of the wall, i.e. e < a. Of course this need not always be the case, 
however, the derivation of the wall equations for e > a is almost identical to e < a, and an 
expansion for general e and a is so notationally cumbersome that we choose to present the 
analysis in this particular case for brevity and clarity. 

First of all, note that the Ricci curvature of the wall is at least of order O(a^), as can 
be seen from (^, since the flat space Ricci curvature is given in terms of products of the 
extrinsic curvature, and the self-gravitating wall has Rab = O(e^). We may therefore set 

Rab = Oi^Rab- (23) 
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Rescaling the extrinsic curvature and parallel coordinates as in ([I2D, and defining 

e 



S 



a 



the equations of motion in the presence of gravity (here characterised by 6) become 



K, = 2aK, 
Kb- 



ab 



(24) 



(25a) 



X" = 2X (X^ - 1 ) - a KX' + a^DaD^X 



K' = -6 



D,K\ - DaK = 26X'DaX 

aR = aKlt, - aK^ + 25 



6{X' - lyhab + a [2KacKMh'"' - Rab - KKab) + 26a' DaXD.X (25b) 

(25c) 

(25d) 

(25e) 
(25f) 



2X'2 + - 1 



-aK\ 



ab 



X'2 - (X^ - if + a^DaXD^X 



To solve these equations, we expand all quantities with respect to a: 



X 

hah 
Kab 
Rah 



Xo + aXi +0^X2 + 



hoab + ahlab + a h2ab + 

Koab + aKiab + Oi' K2ab + 

Rn nb + aRi ah + Ri 



(26) 



iQab -r ani ab ~J~ a 112 ab ~r ' ' ' , 

which also implies similar series for the traces K and R with, for instance, 

K2 = K2abh'^ — Kiabh'f — Koabhf + Kq abhi Jl^i , 



(27) 



iOafe 



where all indices are raised using h^ab- Note that we do not expand in a double series with 
5, since the presence of the a terms in the RHS of ( p5a| ,c,d) means that at any particular 
order in a the series expansion in 5 terminates, as we can see from ( pH]) below. We can now 
solve equations (^) order by order. 
To zeroth order in a we obtain 

(28a) 
(28b) 
(28c) 
(28d) 



(29a) 
(29b) 



Since V{X) is an even function /„ is odd (respectively, even) for n even (respectively, odd). 
As a result Fn is even (respectively, odd) for n even (respectively, odd). By considering the 
derivative of (ESl) we can deduce that 



hoab 


— hoab\o 


Koab 


= Koab\o 


Ko 


= Ko\o- 




= tanhw. 



and (|25t|) is identically satisfied for Xq = tanhw. We define 

Uu)= I duu'^ViXo) 
Fju) = I dufn{u). 
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Rq — Ro\o =^ Roab — Roablo'i (30) 

however, to find the actual value of i?o|o we need to go to first order for X (see equation 
(0) below). 

To first order, we immediately obtain 

hj^ab = h-^ab\o + S-uKpablo — 2SFohgab\o (31) 

and 

- 2X,i3Xf^ - 1) = -X'Jirjo - 35/o]. (32) 
The integrability requirement then constrains Kq\q = 0, and Xi is found to be 

Xi = 35Xo / -— / /oXq = -2 . (33) 

Jo Aq J-oo d cosh u 

Note that Xi is a correction due to the presence of gravity. Then 

i^i = i^ilo - uK^Jo - 45^0^1 - S^Go (34a) 

Klab = Kiab\o — Sfohiablo + SuiTo ac 1 0-^0 ""b 1 — <^ (2/l + Fq) Koab\o 

+ 5Gi{u)hoab- uRoablo, (34b) 
where the functions Gq and Gi are defined by 

Go{u) = y (iw (^4 tanh M sech^wXi + 3/o ) (35a) 
G'i(m) = 2/oFo + y" rfw (4tanhM sech^wXi - 3/o^) , (35b) 

and are both of odd parity. Finally, from ( P5i|) one obtains 

Ro = K'.Jo + 26' [2X',X[ - AXoX.iX', - 1) - Qf^] = K'.Jo - ^5'- (36) 
At second order in a, the equation for X2 is 

VX2 = 6X0X1 - KiX'^ - KoX[ (37) 
= -i^iloX^ + uK'.JoX'o + 6X0X2 + 45X1X^2 ^ ^2^^^/ ^ 

Since all but the first term on the RHS have odd parity, the integrability requirement once 
again constrains Ki\q = 0, and we can solve for X2 giving 

ru ru 

X2 = sech^M / cosh*^-u / {uKlf^\o) secli^udu, (38) 

Jo J-oo 

which is an odd function with respect to u. Therefore, there are no corrections to the Nambu 
equation at first order just as in the case of fiat spacetime. We do anticipate however that 
such corrections will appear at second order, and so proceed to calculate K2. However, as 
we commented earlier, the constraints due to the integrability requirement only pertain to 
the even part of K2, and this in turn depends on the odd part of 
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-25 \2X'^X'^ + + 2X0X2 (Xl - 1) + Xl (3x2 - 1)] (39) 
= 2uKQab\oR'Q'\o + [even terms]. 

Therefore, 

i^2|even = i^2|o + u'K^ab\oKW (40) 

Now, examining the X-equation at order 0{a^) we find 

PX3 = 2X1 + 12X0X1X2 - K^X'q - K^X[ - iTo^s (41) 
and therefore the integrabihty requirement yields the constraint 

-^2|o = — r -^oablo -Rq^Io = - — 1) -^oafe|o-Ro''lo- (42) 
/o(oo) V 6 / 

which gives us the first perturbation of the Nambu equation. 

To summarize, the first corrections to the Nambu equations of motion (reversing the 
rescahngs performed) appear at second order in a and are: 



i^|o = - — -1 Kab\^R%. (43) 




Although we obtained this result assuming e < a, it is in fact quite general, since in a 
similar calculation for e> a rescaling with respect to e gives exactly the same result. Clearly 
from this correction has the correct fiat space limit (p2D, and it would seem that the 
inclusion of gravity simply modifies the second order correction to the equations of motion, 
rather than causing a completely new correction to appear. Indeed the three dimensional 
Ricci curvature can be seen in (|^) to relate to the extrinsic and background geometry of our 
spacetime. 

A totally geodesic {Kab\Q = 0) solution such as the plane-symmetric wall always trivially 
satisfies (^31). What is of real interest, however, is the existence of non-totally geodesic 
solutions verifying the perturbative equations of motion (0). A better understanding of 
the Nambu correction term will involve the computation of the three dimensional curvature 
tensor -Roablo- This task is undertaken in the next section for the specific example of a 
collapsing spherical domain wall. 



IV. THE COLLAPSE OF A SPHERICAL DOMAIN WALL 



In this section we apply the general equations of motion ( ^31) to the case of a collapsing 
spherical domain wall. This is perhaps the simplest non-trivial example of a curved domain 
wall in curved space-time, i.e. where both of our perturbation parameters a and e are not 
zero. The spherical domain wall has been already studied in different contexts using Israel's 
thin wall formalism |[r^ , p!5|j23[| . A first attempt to study the thick case was undertaken in ||20 



but only using equations (P2D valid for a fiat space-time background. 
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Consider, in a spherical system of coordinates {t,r,6,(f)), a non-static scalar field repre- 
senting a domain wall matter coupled to a spherically symmetric metric 



ds^ = A^{t, r)de - B\t, r)dr^ - r^dnjj. 



(44) 



The field equations one has to solve are the coupled Einstein and scalar equations, which 
can be written in a convenient way as 



{ABY 
AB^ 



A" A'B' 



B 



BA 



AB^ 



+ 



B 
B 
2A' 



— = trX'X 



AB^ BA'^ BA^ rAB^ 
DX + 2X(X2 - 1) 



2e^ - eV{X) 
0. 



For small values of e <^ 1, consider the field expansion, 

X = Xo + eXi + 0{e^) 
A=l + eA,+ 0(e2) 
B = l + eBi + 0(e^). 



(45a) 
(45b) 

(45c) 

(45d) 
(45e) 



(46a) 
(46b) 
(46c) 



In order to solve (|45| ) perturbatively one has to first solve the scalar equation ( ^5e| ) to zeroth 
order in e. Then integrating out the Einstein equations ([45a| - [45c|) we can obtain the first 
order e corrections for the metric Ai and Bi. Let us first define the wall's position and 
discuss some general features about the solution before solving (^Sf). 

The location of the wall is defined by the zero of the Higgs field, and will have coordinates 
X" = [t , R{t) , 6 , (p) . We start by making the observation that in order for the wall to be 
identifiably a topological kink solution, R{t) ^ 1, and without loss of generality, we can 
assume that X < (respectively, X > 0) for r < R{t) (r > R(t)). (Note that these are not 
the Gauss-Codazzi gauge coordinates centered on the wall's core.) 

We consider the following initial conditions which are compatible with the fact that, due 
to its spherical symmetry, the wall is not radiating: 



R{t = 0) = R{t = 0) = 0. 



(47) 



The wall's initial bending parameter is defined as a = l/Ri- 

For an inertial observer situated outside the wall, r > R{t), so by Birkhoff's theorem the 
exterior metric is Schwarzschild: 



ds' 



1 



2GM\ , 2 
r J 



1 



2GM\ 



dP + r'^dVL 



2 

Hi 



r > R{t). 



(48) 



In the same way, inside the wall the spacetime metric is flat to a very good approximation. 
The wall's mass, as measured by a distant observer, is M = AnaRf and the Schwarzchild 
radius is given by rs = 2GM. In order for Ri > rs, one must impose the relation a > e. 
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Keeping these considerations in mind, let us now proceed with solving our field equations 
order by order in e in the region close to the wall's core, r ^ R{t). In order to do so, we 
must also expand the function R{t), determining the wall's position, in powers of e. 



R{t) = Ro{t) + eRi{t)+0{e' 



(49) 



Of course Ro(t) will in fact be given in terms of a power series in a: Ro(t) = 7lQ{t) + a^7l2{t) 
etc. 

To zeroth order in e, i.e. in a fiat background spacetime, we can define the unit, exterior 
pointing, normal vector to the wall's core at r = -Ro('^) as 



n 



(i?o,l,0,0)/Vl-i?o • 



(50) 



The intrinsic metric and extrinsic curvature components of the wall (at r = -Ro(^)) are easily 
calculated in turn using (Rf). The extrinsic curvature for e = is then given by 



' (1 - Roy/' 

and the equation of motion (using (|2^)) is 



(51) 



RqvI — Ro 



Ro = ~ Ro 



^-1 

6 



Rn 



[i-R'o 



2\3 



2 



This can be solved iteratively, giving to leading order 



(no 

\Ri 



I -no 



(52) 



(53) 



which can be solved analytically and numerically, as can the 0(a;^) correction (see and 
2^ , |25[| ), giving the wall trajectory as shown in figure |I]. To 0{a') we note that Rq satisfies 



1-^0-;^ 



,2C Rr 
Rf V Rl 



(54) 



we note that 4^Ro = 0(a" ^] 



where C = /2(oo)//o(oo) = — 1. From ( |^ ) and 

In order to solve for the spacetime metric, we need Xq. From the work of the previous 
section we know that Xq = tanhw + X2, where X2 is given from equation (120) as 



X2 



Rli 



6 

r— sech^u / cosh'^ u u sech'^u, 

l-Rl) Jo J-oo 



and is clearly of order O(a^). Here, u is the proper distance from the wall. Using 
note that 



t = t' 



r = Roit') 



uRo{t') 



1 - Ro it') 
u 



1 - Ro it') 



(55) 

we 

(56a) 
(56b) 
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FIG. 1. The trajectory for the collapsing wall in flat spacetime. The solid line is the leading 
order Nambu trajectory, and the dashed line the trajectory including O(a^) corrections. For 
illustrative purposes the rather large value of a = 1/5 has been chosen. 



where (t', u) are the coordinates of the point [t, r) in Gaussian Normal gauge. 

In order to find the metric, first note that we can directly integrate ( [45 a|) and ( [45b|) 
implicitly, finding 



Bi = — 
2r 



1-Rl{t') 



2r Jo 



+ viXo) 
] X',' + V{Xo) 



dr 



dr + 



l + Rl{t') 



X'^dr, 



(57a) 
(57b) 



where we have used (H) to obtain u^ + u'"^ = [1 + Rl{t')]/[1- Rl{t')]. We may now substitute 
the form of Xq{u, t') to the required order, replace the r-integral by a u- integral along a line 
t = const, then Taylor expand t' around m = to the required order in a. For example, Bi 
gives the mass function via = 1 — 2GM{r) /r =^ GM{r) = erBi. Computing Bi from 
the above expression to O(a^) yields 



GM{r) 



^&^h{u) + 2e^J,{u) + 6i?f sech\x2(«) + ef^ \h{u) + QRluUn)] (58) 



where we have put 



fn{u) 



c/mm"\/(Xo). 



(59) 



This gives the ADM mass (using ^4|) as 
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10 15 20 25 



FIG. 2. The r — r component of the metric for the collapsing wall solution of figure |l] for various 
values of t up to tc, the critical value of t for which the wall reaches its Schwarzschild radius. Again, 
for illustrative purposes, the rather large value of e = 1/14 has been chosen. 

~ ~ 4 / TT^ \ 

GMadm = lim erB, = -—±=f,{oo) + 6-^/2(00) = -e /^f + - - 1 , (60) 



which is indeed constant, and agrees to leading order with the expected result. Figure |^ 
shows the evolution of the grr metric component up to the formation of a black hole. 
The gravitational Nambu equation for e 7^ at u = is given by 

R + ^nA' + {A' A + RBB)n + (f " ^) + ^ (| - ^) = 0(«'e), (61) 

where n = B~^ — R^A~^. Note that the gravitational correction to A and B is 0(5) not 
0(e), where 5 = e/a was defined in (^3p. This means that gravitational corrections to the 
wall motion potentially could appear at 0(5). Since the flat space wall equations are intially 
0(a), we could have the catastrophic situation that adding gravity swamps the wall motion, 
and has a superdominant effect. In fact, this turns out not to be the case. The computation 
of the A and B contributions to 0(5) in ( |6lD shows that they cancel. Instead, the leading 
order correction appears at 0(5a; = e), which is therefore subdominant to the flat space 
motion, and is plotted in figure ^ 

R = -1(1 - R^) + 2e(l - i?2)3/2 + 26(41n2 -l)j^^ ^^ _ ^ ^g^) 
R 3R 

Note that we have included the 0{ae) correction, since it is the same order as the finite 
width correction, however, this has not been used in the computation of the corrected 
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Nambu trajectory in figure ^. We should stress that this correction does not mean that the 
wall trajectory is no longer Nambu, it simply gives the right gravitational alteration to the 
trajectory to allow the wall to remain a minimal surface in the now curved spacetime. What 
is important however, is that the correction is subdominant to the flat space motion and is a 
simple 'nudge'; there is no evidence that adding gravity gives any catastrophic effect which 
could force a trajectory to be totally geodesic. 




1 2 3 4 5 

t 



FIG. 3. A comparison of the leading order (i.e. Nambu) flat space bubble trajectory (solid line) 
and the leading order (Nambu) trajectory in the presence of gravity. As before, for illustrative 
purposes the rather large values of a = 1/5 and e = 1/14 liave been cliosen. Although the bubble 
actually reaches its Schwarzschild radius at t ~ 2.45 (see fig |2|), we have chosen the i-range for 
direct comparison with figure |l|. 

What is of real interest to us however, at least in the context of our general equation 
(^), is the correction to the Nambu motion. Indeed, having a particular example, we can 
explicitly calculate Rat from Gauss's equation (|^). Since we have not rescaled our quantities 
here, we can check that the correction term is indeed of the right order, and get an idea 
of the physical implication that the correction induces (at least in this particular example). 
We obtain. 

The first term is a correction term of order 0(a'^) due to the bending of the wall, already 
present in flat space-time as predicted in [Q. The second correction term is due to self- 
gravity appearing at 0(ae), as we were indeed expecting (for a spherical wall a > e). Note 
that the finite width gravitational correction can dominate the extrinsic curvature correction, 
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(63) 



depending on the relative magnitudes of e and o? . Finally from the overall positive sign 
we can deduce that the corrections to the Nambu motion induce a slowdown of the wall's 
collapse. 

Using (|62|) and (|63|) we can write the full equation of motion (i.e. including the dominant 
deviation to Nambu motion) for a collapsing spherical wall up to and including order ae, 

R= - ^(1 - R^) + 2e(l - Ri^fl^ + ^ — — - — (1 - 6i?2 + 3i?^) 

R oR 

+ (^ -6)o-^— rf- (64) 



Before leaving the collapsing wall, it is worthwhile comparing our thick wall calculation 
with the results of Ipser and Sikivie, [|1^], obtained for the collapsing thin wall. In order 
to obtain the Israel thin wall approximation one should reintroduce the width parameter w 
which has been set equal to unity throughout this analysis. Then a = w/R (where R is a 
typical radius of curvature of the wall) is a dimensionless parameter. The Israel limit is then 
obtained by letting a and e tend to zero, while keeping their quotient 6 fixed. This amounts 
to keeping the horizon distance from the wall finite i.e., keeping the wall self-gravitating (as 
one should in this formalism). This also gives us the correct limit without involving the 
normal coordinate u in the limiting procedure. Taking this limit in (|64[) gives 



i? = -4(1 -i?') + 2(5(1 -i?2)3/2 (65) 
R 

which is obviously only correct to 0(5). The translation to the calculation of Ipser and Sikivie 
is not direct, since our coordinates correspond to the interior coordinates of their bubble, 
however, once the correct correspondence is made, we do indeed find precise agreement to 
order 0{5)=0{Gcr). 



V. CONCLUSIONS 

We have obtained a general equation of motion (^31) describing a moving wall in curved 
spacetime. This was achieved by analytically solving the Einstein and scalar matter field 
equations order by order with respect to two parameters: the wall's bending parameter, a, 
and the gravitational strength parameter, e, expressing the curving of spacetime. We then 
considered a particular example, the collapsing spherical domain wall, which is perhaps the 
simplest non-trivial example with both parameters different from zero. In the context of 
this example we found that the corrections to the flat space Nambu motion tended to slow 
down the wall's collapse. 

Throughout this paper we have considered a thick wall with scalar (Higgs) matter. This 
was done in order to examine the problem in its analytic (with respect to the spacetime 
metric) and most general context and also in order to pick up finite width gravitational 
corrections. As we discussed at the end of the previous section, in order to obtain the Israel 
thin wall approximation one takes the limit a — * and e — > keeping 6 fixed. Although the 
spherical wall was only explored to 0(5), the arguments in section |T| show in all generality 
that the corrections to the Nambu equation of motion K = are finite width corrections. 
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and hence disappear in the thin wall limit, giving simply gravitational corrections to the flat 
space trajectory. 

A totally geodesic {Kab\o = 0) solution such as the plane-symmetric wall always trivially 
satisfies (|43|). With respect to the physical motion of the wall a totally geodesic solution 
is trivial and in particular a totally geodesic wall will not emit gravitational waves. What 
is of real interest physically, is the existence of non-totally geodesic solutions such as the 
spherical domain wall satisfying (^). There have been claims however (see for the related 
example of a cosmic string, p6[, and references therein), that the presence of gravity in 



general constrains a defect's core to be totally geodesic. For a domain wall in particular it 
has been claimed that the presence of gravity induces the wall to lose its dynamical degree 



of freedom and not to radiate We should stress that throughout our treatment this 

constraint has not appeared and gravity affects the motion of the defect in a very natural way 
in the sense that corrections to the Nambu motion appear at the same order as in the case of 
a flat background spacetime. We suspect that the reason for this discrepancy is that in |^ 
rather speciflc asymptotic boundary conditions have been placed on the spacetime, namely 
that it asymptote the static planar domain wall solution. If, however, a wall is oscillating, we 
expect that its effective equation of state will change, analogous to that of the wiggly cosmic 



string |2^, which will change the spacetime metric even asymptotically. Mathematically, 
this can be seen via the divergent behaviour of the metric perturbations in the Gaussian 
Normal gauge due to proper motion of the wall, and was discussed in the context of higher 



dimensional domain walls in |29]. Following [2q|, we expect that the equation of state of 



the wall will have the form cxT^ = ctq, or, that perturbatively the effect on the energy- 
momentum tensor of the wall will be to increase the energy by Sa, and decrease the tension 
by Sa/2. Note that this perturbation is tracefree, and localised on the wall. In fact if we 
regard our domain wall as a 2 + 1-dimensional universe, this energy momentum would be 
that appropriate to a radiation dominated cosmology. Such systems have been explored in 
a different context in, for example, [^. 

As a flnal point it should be clear that although we have considered a four- dimensional 
spacetime our equations of motion are valid for any [n — 2)-brane of an n-dimensional 
spacetime. Our analysis however relies heavily on the fact that the wall is a hypersurface 
of the ambient spacetime i.e. that codimension is one. If we were to consider dynamics of 
strings for instance, the picture could in principle be quite different. 
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